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ABSTRACT 


In this paper the data induced from lightning is synthesized using suitable inversion techniques so as to delineate 


imperative information regarding availability of groundwater. 
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INTRODUCTION 


Earth profiling is an important task being undertaken all over the world for delineating feasible sites of minerals, oil or 
groundwater occurrence. Experts and amateurs working in this field are practicing numerous methods for exploring sub- 
surface information. Delineating mineral resources, groundwater and the reduction of the risk posed by geological hazards 
require more and more accurate and detailed knowledge and mapping of the underground (Fancsik, T. et. al., 2021). Here, 
in order to extract relevant information regarding the occurrence of groundwater, lightning data is worked upon on the lines 
suggested by Grimm, 2002. The techniques of generalized linear inversion and Marquardt inversion are discussed with the 
nature of their applicability (Evans et. al., 2002, Haber, 2004) in context of the problem of earth profiling. The essence of 
inversion technique is that the processing and interpretation of data of complex subsurface geological structures is 
performed by discretization and the inversion procedure formulated for expansion coefficients. The method is found 
suitable for groundwater exploration as even with the introduction of relatively few expansion coefficients, a suitable 
resolution can be achieved so that the problem to be solved leads to an over determined inverse problem. Also, the inverse 
and forward modelling leads to delineating not only sites for groundwater occurrence but also thickness of sub-surface 


layers (Menke, 1984; Dobroka et. al., 2016) 
THEORY 


The fundamental principle of current propagation forms the basis of this study (Diendorfer, 2015; Ushio, et al. 2015; Li, D. 
et al. 2016). The lightning current induced in the earth provides valuable information regarding earth sub-surface profile. 
The data is incorporated to develop a forward problem to see through the sub-surface (Dweight, 1936; Goldman et. al., 
2004). Establishment of a forward problem requires resistivity kernel function to be approximated as an exponential 


function of unknown coefficients, Koefoed, 1976. 


K)= Shee 
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The mathematical elaboration of this kind expresses the variation of resistivity of different layers beneath the earth with 
respect to their extension in the downward direction. Here, f; is the coefficient of exponential function; €; is the exponent 


parameter of kernel function and j is the number of point of observations on the earth surface. In matrix notation 


[Kx = [E] mxplFlpx1 (2) 

e7e4s : e7Fp4s 

e ish 
[2] = 
—,Ay-: 
i= au a : 
ee Fs4m : 2 : e *p4m (3) 
The unknown coefficient matrix can be obtained as follows 

[F] x1 = (G] pxm [U] ux: (4) 


Here G is the generalized linear inversion (Backus et al., 1967; Snider, 1991) of function G which delimits the 


effect of matrix being singular and non-square dimensions. 
(5) 


The objective behind employing generalized linear inversion is that, it offers a delicate solution of the non-linear 
system problem (Backus et al., 1967). Formulating a detailed procedure for obtaining inverse of a non-square matrix, the 


GLI is obtained as 
[G] = [G7 GyG? (6) 


In order to assess the applicability of the method a three layer model is assumed with the study area extending up 
to a range of 5 Km. In order to simplify the complexity of problem the study area is decomposed into a set of 12 
observation points. Table | beneath provides computed values for equation (6) using equation (5). Coefficients of kernel 


function, and then finally take on values as given in Table 2. 


Figure 1 show once the exponential function and its coefficients are estimated from the lightning induced 
potential difference data computed on the surface of earth, the next procedure starts with the evaluation of kernel function 


for different observation points. 


Matlab tool is utilized to get a discreet view of the variation of Kernel function as obtained through computed 
values. It is observed that the resistivity kernel function incorporates the characteristic of in homogeneity on account of 


layer structure and acquires unique values for different layers and different observation points. 
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-741923 | 5851848 | -38553470 | 26428766 | 25950092.37 | 16389112.96 
904672.8 | -7436498 | 48993629 | -33585717 | -32977347.64 | -20827254.5 
-2043136 | 1670835 | -11009428 | 7547487.2 | 7410643.959 | 4680253.824 
10463..96 | -86884.64 | 574041.37 | -393845.1 | -386628.991 | -244156.435 
-78.8455 | 705.1738 | -4819.072 | 3344.567 | 3274.10838 | 2064.81341 
69714 | -6.196523 | 4691145 | -35.883 -34.29348 -21.32920 
7740066.09 | 1128069 | -3785906 | -7457099 | -10240706 | -23096124 
-9836062.6 | -1433554 | 4811109 | 9476449 13013847 29350461 
2210321.103 | 322122.9 | -1081165 | -2129549 | -2924464 -6595581 
-115288.156 | -16775.2 | 56437.37 | 1111133.3 | 152605.11 | 344129.94 
972.627 138.1164 | -481.9275 | -945.089 -1298.24 -2917.657 
-9.7681673 | -.97143 | 5.551439 | 10.4183 14.10544 31.106661 
Table 2 
I] F@ 
1 | .2857 
2 | 164642 
3 | 14285 
4 | .08231 
5 | .07142 
6 | 041155 


Figure 1: Resistivity Kernel Function for Observation Points. 


Table 3 Shows the resistivity kernel function as estimated is utilized for the interpretation of resistivity transform 


function using a recurrence relation based on the work of Pekeris, 1940 and later extended by Koefoed, 1976. 


Ty Q) = py [2K Aj) +1] (7) 
Table: 3: Resistivity Transform Obtained Using Lightning Data 
AQ) 1 5 333 25 2 .1666 14285 
T(A) | 5.860585 | 5.860599 | 5.860975 | 9.862659 | 9.866288 | 11.871897 | 11.879236 
Table 4 
125 111 ll 05 0333 
11.888007 | 13.897951 | 13.908863 | 14.047349 | 14.204395 


Simultaneously, an artificial three layer earth model is assumed with parameters 
p,=100;h = 10 


p14 =4;h= 30 


www.iaset.us editor @iaset.us 


32 Dr. Jyotika 


For this assumed model too resistivity transform is evaluated using the recurrence relation given by Pekeris, 1940. 


“po T w(A)t+Ty-1i(A) 
Tyo A) = ra ty (8) 
so 
i Pn-— 4/2 
Where, 
T’ n(A) = Pn (tanh A hy) (9) 


Table 5: Resistivity Transform 
Obtained for Assumed Model 


a Tmo) 

1 9.411091 
Pe) 11.93513 
333 13.31036 
25 14.22263 
2 14.894 
.1666 15.8404 
14285 15.84604 
125 16.20277 
111 16.50719 
1 16.77146 
OS 18.35025 
.0333 19.1764 


Observed and Computed Resistivity Transform for different points of observation 
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Figure 2: Resistivity Transformation. 


Figure 2 shows the resistivity transforms functions as depicted above are a non-linear function of model 
parameters, p and h. A discreet and iterative inversion scheme is sought out for the complete resolution of earth data which 
is further refined to get closer to groundwater locations. Koefoed (1976) proposed an interactive iteration technique for 


bridging the gap between the assumed earth model and the computed parameter values using the formulation 
wae = . eTi 
Ty =Tyot L(P1— Pua) 55, (10) 


oTi 


Ty — Tuo = XP, — Puro) Ty 


tL) 
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With, 


being called as the sensitivity matrix. 
Also, 


(P; — Puyo) = AP 

Is considered as the correction key for each iteration step 

Then the equation in matrix form would become 

[AT]=[4P] JU] (12) 
This correction step is computed while inverting the sensitivity matrix as 

[AP] = [/] [47] (13) 
Here j is the generalized linear inversion obtained as 

r=? (14) 
Then for each iterative step a new predicted data set will be 

T new = Tuo(A)+ OP (15) 
The steps goes on this way until the assumed parameters approach the true data and thus with the final parameters 


is the accepted one corresponding to the observed data set. 


The above discussed generalized inversion of sensitivity matrix and the parameters F and G results as a 
consequence of direct attempt of minimizing cumulative square of data matrix (Backus et al., 1967: Turai, E., et al., 2001; 


Tavakoli, S. et al., 2016; Gyulai, A., et al., 2017). An acute problem arises when the inverse of | does not exist. In 


such cases, an alternate route is approached by Levenberg (1944) and later by Marquardt (1963). The procedure named as 


the damped least square inversion is later discussed in detail by Inman, 1975. 


An extensive analysis of both inversion techniques reveals that there is absolutely no difficulty with the use of ridge 
regression or the least square solution. However, it is observed by various authors that a perfect resolution of the assumed model 


is not obtained. Thus a generalized linear inversion technique has been used for the detailed analysis of the problem. 
CONCLUSIONS 


The present paper discusses an efficient method for the prediction of earth resistivity data utilizing the information 
provided by lightning phenomenon. The computation procedure is simple and straight forward with the added significance 
that the lightning data is directly used in the matrix method. The inversion techniques discussed here rely on the condition 
of whether the problem is underdetermined or over determined. MS-Excel 2007 and Matlab are customized for fast 


convergence of results. 
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